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ABSTRACT 

We present time-resolved optical spectroscopy of the afterglows of the gamma-ray bursts GRB 990510 
and GRB 990712. Through the identification of several absorption lines in the first epoch GRB 990510 
spectrum, we determine the redshift for this burst at z > 1.619. No clear emission lines are detected. 
The strength of the Mgi feature is indicative of a dense environment, most likely the host galaxy of 
GRB 990510. Although the host is extremely faint (V ^ 28), the GRB afterglow allows us to probe its 
interstellar medium and - in principle - to measure its metallicity. The optical spectrum of GRB 990712 
(whose host galaxy is the brightest of the known GRB hosts at cosmological redshifts), shows clear 
features both in emission and absorption, at a redshift of z = 0.4331 ± 0.0004. On the basis of several 
line emission diagnostic diagrams, we conclude that the host galaxy of GRB 990712 is most likely an 
Hll galaxy. We derive a unreddened [On] star formation rate of 2.7 ± 0.8 M© yr~^. Correcting for 
the measured extinction intrinsic to the host galaxy (Ay = 3.4ti'7), this value increases to 35I25 ^o 
yr~^. The [Oil] equivalent width, compared to that of field galaxies at z < 1, also suggests that the 
host of GRB 990712 is vigorously forming stars. We employ the oxygen and H/3 emission-line intensities 
to estimate the global oxygen abundance for the host of GRB 990712: log(0/H) = -3.7 ± 0.4, which is 
slightly below the lowest metallicity one finds in nearby spiral galaxies. For both GRBs we study the 
time evolution of the absorption lines, whose equivalent width might be expected to change with time if 
the burst resides in a dense compact medium. We find no evidence for a significant change in the Mgii 
width. 
Subject headings: gamma rays: bursts — galaxies: distances and redshifts — cosmology: observations 



1. INTRODUCTION 

In February 1997, the Italian-Dutch satellite BeppoSAX 
enabled a breakthrough in the understanding of gamma- 
ray bursts (GRBs) by providing an accurate position for 
the prompt X-ray emission of a GRB. This led to the dis- 
covery of the first X-ray afterglow of a GRB (Costa et al. 
1997), and, independently, to the identification of the first 
optical counterpart of a burster (Van Paradijs et al. 1997). 
Since then, several X-ray, optical and radio counterparts 
of GRBs have been detected. These multi-wavelength af- 
terglow observations can be explained reasonably well by 



simple fireball models (for recent reviews see Piran 1999, 
and Van Paradijs, Kouveliotou & Wijers 2000). GRB dis- 
tance determinations are crucial in the effort to establish 
the physical nature of their progenitor (s). The observed 
redshift distribution of the 'normal' afterglows (i.e., ex- 
cluding GRB 980425, which is associated with supernova 
SN1998bw at z = 0.0085; Galama et al. 1998), ranges from 
z = 0.43 (Galama et al. 1999, and this paper) to z = 3.42 
(Kulkarni et al. 1998). 

Although major advances in the understanding of GRBs 
have been made over the past few years (thanks to the 
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detection of afterglows), the physical nature of their pro- 
genitor(s) remains unclear. The most popular models are 
(i) the collapse of a rotating massive star (Woosley 1993, 
MacFadyen & Woosley 1999) and (ii) the merging of two 
neutron stars, or a neutron star and a black hole (Narayan, 
Paczyhski & Piran 1992, Janka et al. 1999). The former 
('collapsar') model has trouble producing GRBs with du- 
rations shorter than a couple of seconds and predicts that 
every afterglow is accompanied by a supernova of a type 
(Ic) similar to SN1998bw (MacFadyen & Woosley 1999). 
Furthermore, in the collapsar environment it is likely that 
the optical light of the afterglow is heavily absorbed by the 
surrounding dusty medium. Given the massive progeni- 
tors, one expects the frequency of GRBs to be strongly 
correlated with the cosmic star formation rate; the latter 
remains one of the great unresolved issues in astronomy of 
today. The compact star merger scenario can make short 
GRBs as well as long ones (although a 10^^ G magnetic 
field is probably needed for the latter, see Meszaros 2000); 
some of these mergers are expected to occur in low-density 
environments, possibly located several kiloparsec outside 
their host galaxies. This is due to the large kick velocities 
imparted to the compact objects from the two respective 
supernovae (~ 250 - 300 km s~^; Hansen & Phinney 1997), 
combined with the long time between the birth of the sys- 
tem and the merger occurrence (10^ - 10^ years; Portegies 
Zwart & Yungelson 1998). Consequently, since the opti- 
cal afterglow brightness depends on the density of the cir- 
cumsource medium, some of these bursts may not show an 
afterglow at all. These could account for the 'dark' burst 
population, i.e. bursts for which only X-ray afterglows 
have been found. One way to discriminate between these 
two models is by studying the immediate environment of 
the burst. In the collapsar model the circumsource den- 
sity is expected to drop with distance as r~^, due to the 
expanding stellar wind of the SN progenitor, while in the 
binary merger scenario a constant, relatively low-density 
ambient medium is most plausible. 

If the GRB source resides in a compact, gas-rich environ- 
ment (which is expected in the collapsar scenario), the af- 
terglow spectrum might show time-dependent absorption 
features (such as Lya and Mgii) due to the gradual ion- 
ization of the surrounding medium (Perna & Loeb 1998). 
In this case a decrease of the absorption-line equivalent 
widths (EWs) with time is expected. On the other hand, 
spectroscopic observations of the star HD 72089, situated 
behind the Vela supernova remnant, show an increase of 
an order of magnitude of the absorption strengths of ele- 
ments such as Al and Fe, over the velocity range spanned 
by absorption in the remnant (Jenkins & Wallerstein 1995; 
see also Savage & Sembach 1996). This is attributed to the 
destruction of the dust grains, due to the propagation of 
the SN shock, which causes the release of elements (such 
as Fe and Mg) that are frozen in the dust. Thus in a dusty 
environment that is being 'shocked' by a GRB explosion, 
one might expect the strength of the absorption lines to 
increase in time. 

In order to test these theories of GRB genesis, and 
the effects of GRBs on their environments, we have an 
on-going program to obtain spectra of GRB afterglows 
using the Very Large Telescope (VLT) of the European 
Southern Observatory (ESO) at Paranal, Chile. Here we 
present results on two of these bursts: GRB 990510 and 



GRB 990712. 



1.1. GRB 990510 



GRB 990510 was observed on 1999 May 10.36743 UT 
with the Wide Field Camera (WFC) unit 2 onboard Bep- 
poSAX, which localized the burst at R.A. = 13''38"06", 
Deck = -80°29'30" (J2000.0), with an error radius of 3' 
(Dadina et al. 1999). The BeppoSAX Gamma Ray Burst 
Monitor (GRBM) recorded an 80s event with a multi-peak 
structure (40-700 keV) . The average and peak intensity in 
the WFC unit 2 (2-28 keV) was about 0.7 and 4.3 Crab, 
respectively. The burst position as determined with Bep- 
poSAX is consistent with that of the Inter- Planetary Net- 
work (IPN; Hurley et al. 1999), using the Burst And Tran- 
sient Source Experiment (BATSE) onboard the Gompton 
Gamma Ray Observatory, and Ulysses. BATSE recorded 
a fluence above 20 keV of 2.56 x 10"^ erg cm^^ (Kippen 
et al. 1999), ranking it in the top 9% of the burst fluence 
distribution. 

We discovered the optical afterglow on images taken at 
the South African Astronomical Observatory (SAAO) Im 
telescope (Vrecswijk ct al. 1999) and subsequently trig- 
gered our VLT program to take spectra and polarimetric 
images. Here, we present the time-resolved spectroscopy. 
Our polarimetric observations resulted in the first polar- 
ization detection of a GRB afterglow (Wijers et al. 1999; 
see also Covino et al. 1999), while our photometric obser- 
vations show an achromatic break in the BVRIJHK light 
curves, which is most likely due to the burst emission be- 
ing coUimated (Rol et al. 2000a; see also Stanek et al. 
1999, and Harrison et al. 1999). Fruchter et al. (1999) 
have used HST to estimate Vhost ^ 28, and do not find evi- 
dence for a supernova (SN) of the same type and brightness 
as SN 1998bw in the late-time hght curve of GRB 990510. 
Recent HST observations (April 2000) appear to detect 
a faint galaxy (V ~ 28) at the position of the early op- 
tical transient, which, if real, is most likely the host of 
GRB 990510 (Fruchter et al. 2000, Bloom et al. 2000). 

1.2. GRB 990712 

GRB 990712 triggered the GRBM and WFC unit 2 on- 
board BeppoSAX on 1999 July 12.69655 UT. The burst 
lasted for about 30s, had a double-peaked structure, was 
moderate in 7 rays, and was accompanied by one of the 
strongest prompt X-ray counterparts ever observed (Heise 
et al. 1999). The WFC unit 2 located the burst at R.A. 
= 22''31'"50", Deck = -73''24'24" (J2000.0), with an er- 
ror radius of 2'. Unfortunately, neither flux nor fluence 
levels are reported in the literature. Again the SAAO Im 
telescope was successful in hunting down the GRB after- 
glow (Bakos et al. 1999), which allowed us to quickly alert 
the VLT staff for spectroscopic, polarimetric and further 
photometric follow-up observations. The host galaxy of 
GRB 990712 is the brightest of the known GRB host galax- 
ies, with R = 21.8 and V = 22.3 (Sahu et al. 2000). The 
VLT polarimetric images exhibit a signiflcant degree of 
polarization of the afterglow of GRB 990712 which seems 
to vary with time, while the polarization angle does not 
change with time (Rol et al. 2000b). These observations 
cannot be easily reconciled with afterglow polarization the- 
ories. The photometric measurements show a common 
temporal power-law decay of the transient source, over- 
taken by the bright host galaxy at late times; no evidence 



Table 1 
Log of the GRB 990510 and GRB 990712 spectroscopic observations. 



UT date (1999) 


time since burst 
(days) 


grism 


wavelength range 

(A) 


resolving power 


S/N at 6500 A 


GRB 990510 












May 11.179 


0.811 


G150I 


3700-7700 


185 


91 


May 11.203 


0.836 


G150I+OG590 


6000-9000 


280 


78 


May 12.123 


1.755 


G150I 


3700-7700 


185 


6 


May 12.146 


1.779 


G150I+OG590 


6000-9000 


280 


6 


May 14.273 


3.906 


G300V and G300I 


3880-9255 


420 and 680 


9 


May 16.254 


5.887 


G150I 


3700-7700 


185 


5 


GRB 990712 












July 13.182 


0.485 


G150I 


3700-7700 


185 


26 


July 13.421 


0.725 


G150I 


3700-7700 


185 


21 


July 14.181 


1.485 


G150I 


3700-7700 


185 


14 



Note. — The total exposure time for all spectra was 1800 sec. 



is found for a supernova of type SN1998bw (Sahu et al. 
2000; see also Hjorth et al. 2000). 

The organization of this paper is as follows: in §2 we 
present the observations and data reduction methods. In 
§3 we display and discuss the spectra of GRB 990510, fol- 
lowed by GRB 990712 in §4. We study the absorption-hne 
intensity evolution in time for both bursts in §5 and de- 
scribe our conclusions in §6. 

2. OBSERVATIONS 

After the optical identification of both GRB 990510 
(Vreeswijk et al. 1999) and GRB 990712 (Bakos et al. 
1999), we triggered our VLT target-of-opportunity ob- 
servation program and obtained several low-resolution 
spectra at various epochs with the FOcal Reducer and 
low-dispersion Spectrograph (FORS), mounted at the 
Cassegrain focus of the ESO VLT-UTl Antu telescope. 
The date of observation, grism used, wavelength range, re- 
solving power and the signal-to-noise at 6500A are listed 
in Table m A slit width of 1" was used for all spectra. 
Grism G150I approximately covers the wavelength range 
3700-9000 A. However, redward of 6500 A the second order 
starts to contaminate the first order (cf. FORS User Man- 
ual 1.3). To obtain a clean spectrum longward of 6500 A 
for GRB 990510, we also took spectra with an order sep- 
aration filter (OG590), using the same grism. However, 
due to the low sensitivity of the CCD shortward of 3700 
A, the impact of the overlap is negligible shortward of 7700 
A. Therefore, we have summed the blue and red spectra 
over the region 6200-7700 A and combined this part with 
the single blue and red spectra into a continuous spectrum 
over the entire wavelength range. For the grisms G300V 
and G300I, we have simply connected the blue and red 
parts into one spectrum. For GRB 990712, we used grism 
G150I without order separation filter, and thus these spec- 
tra are usable over the wavelength range 3700-7700 A. 

The raw spectra were bias-subtracted, and flat-fielded 
with a normalized combined set of lamp flat-fields. Sub- 
sequently, cosmic rays were removed interactively along 
the afterglow spectrum and each sequence of images was 
summed into combined images. The spectra were opti- 



mally extracted from the combined images, and wave- 
length calibrated using a standard Helium-Neon-Argon 
lamp. The r.m.s. error in the wavelength calibration is 
approximately 0.25 A. Since no useful spectra of standard 
stars were taken neither during the first two nights for 
GRB 990510 nor for GRB 990712, we have flux-cahbrated 
these spectra using the BVRI light curve data of Rol et al. 
(2000a) for GRB 990510 and the VRI data of Sahu et al. 
(2000) for GRB 990712. The spectra of the other nights 
were flux-calibrated with a spectrophotometric standard 
star, which resulted in flux levels that are consistent with 
the photometry at the same epochs. 

To flux-calibrate the spectrum of GRB 990510, we fitted 
the light curves with a smoothly connected broken-power- 
law model (Harrison et al. 1999, Stanek et al. 1999, 
Beuermann et al. 1999), while for GRB 990712 we used 
a simple power-law model with a host galaxy contribu- 
tion. We determined the magnitudes at the times when 
the spectra were taken (see Table |l|), using the fits to the 
light curves. These values were then corrected for the esti- 
mated Galactic foreground extinction: E(B-V) = 0.2 and 
E(B-V) = 0.03 (Schlegel, Finkbeiner & Davis 1998) for 
the May and July burst, respectively. The E(B-V) value 
is translated into an extinction at a given wavelength using 
the standard Galactic extinction curve of Cardelli, Clay- 
ton & Mathis (1989). The magnitudes were transformed 
to fiuxes (Fukugita, Shimasaku & Ichikawa 1995) and were 
fitted with a power-law spectrum F oc i^''. For GRB 990510 
we find /3 = -0.6 ± 0.1 (first night), and (3 = -0.7 ± 0.1 
(second night). For the three spectra of GRB 990712 we 
obtain /3(1) = -1.1 ± 0.2, /3(2) = -0.9 ± 0.2 and /3(3) 
-0.9 ± 0.2. We have not taken into account the extinc- 
tion intrinsic to the host galaxy; see §4. These slopes are 
quite usual for GRB afterglows. We also fitted the global 
profile (excluding the absorption and emission lines) of all 
the wavelength-calibrated spectra with a 4th-order Cheby- 
chev polynomial. To obtain the flux-calibrated spectra, we 
multiplied the wavelength-calibrated spectra by the ratio 
between the power-law fit based on the photometry and 
the global spectral profile fit. We estimate the error in the 
flux calibration for all nights to be about 15%. 



The equivalent width (EW) of the spectral lines was 
measured using the splot routine in IRAF; we used both 
a Gaussian fit, and also simply summed the difference be- 
tween the pixel value and the continuum for each pixel 
over the line. Both methods gave similar results. The er- 
ror is mostly dominated by the uncertainty in the contin- 
uum level, which was estimated by placing the continuum 
at a high and low level (roughly corresponding to the ± 
r.m.s. value of the continuum in the vicinity of the line). 
We take the mean value of these two estimates as the EW, 
and half their difference as the error. We also calculate the 
Poisson error from the object and sky spectrum, which we 
quadratically sum with the measurement error to obtain 
the total error. We then corrected the EW for the host 
galaxy contribution in case of absorption features (i.e. di- 
vided the measured EW by the afterglow fraction of the 
total light) and for the afterglow contribution in case of an 
emission line (i.e. divided by the galaxy-light fraction of 
the total). The host galaxy fraction of the total light for 
GRB 990510 is neghgible (see §3), while for GRB 990712 
we obtain 0.23 ± 0.03, 0.30 ± 0.04 and 0.47 ± 0.06 for 
epoch one, two and three, respectively, from the V and 
R band light curve fits of Sahu et al. (2000). The error 
was estimated from the different values for the host galaxy 
magnitude as given by Sahu et al. (2000) and Hjorth et al. 
(2000). These errors are propagated along with the errors 
in the EW measurements. Finally the EW is converted 
to the rest frame by dividing by (1-l-z). The emission-line 
fluxes are independent of the afterglow contribution to the 
total light. 

3. THE ABSORPTION-LINE SPECTRUM OF GRB 990510 

We have identified several absorption lines in the first 
epoch spectrum of GRB 990510: Mgii A2800, Mgi A2853, 
Fell AA2344,2383,2600, and possibly Ahii A1683 and 
Crii/Znii A2062. The observed wavelengths, identifica- 
tions, corresponding redshifts, and equivalent widths (in 
the absorber's rest frame) are listed in Table 0. The lines 
and the telluric absorption features are also indicated in 
Figure 1. 

These line identifications can be verified by taking into 
account the oscillator strength, ionization potential and 
relative cosmic abundance. The observed EWs of Fell at 
2344 A, 2383 A, and 2600 A are in reasonable agreement 
with their relative oscillator strengths (Morton 1991). 
Other Fen hues, such as Feii A1608 and Feii A2587 were 
not observed, but these have smaller oscillator strengths, 
consistent with their non-detection. However, the oscilla- 
tor strength of the undetected Aim A1855 is twice that 
of the detected Aim A1863 line, and so should have been 
detected at an observed wavelength of 4858 A. This makes 
the identification of the line at 4883 A with Aim question- 
able. A similar argument can be made for the reality of the 
Crii/Znii A2062 line. Both lines are expected to be accom- 
panied by a stronger partner (Cm A2056 and Znii A2026), 
which is not detected. The presence of Fen suggests that 
the medium probed with our line of sight is most likely 
one of low-ionization, and therefore Aln A1670 is expected 
to be present as well given the similarity between the Fe 
and Al ionization potentials, but again nothing is detected. 
However, we clearly detect Mgii, which is a blend of Mgii 
A2796 and Mgii A2804, and observe Mgi A2853 in absorp- 
tion as well. 



^ erg cm ^ (Kippen et al. 



The weighted mean redshift of the identified lines is 
z — 1.6187± 0.0015. On the basis of this redshift we can 
possibly identify Feli(2250) at 5893 A, although we do not 
expect to detect this line due to its low oscillator strength. 
Since no clear emission lines arc detected, this redshift is 
certainly a lower limit to the redshift of the GRB afterglow. 
However, the strength of Mgi A2853 suggests we are prob- 
ing a dense, low-ionization medium, very likely that of the 
host galaxy. Assuming isotropic emission, z = 1.619, Hq = 
70 km s~^ Mpc~^, rio=0.3, and A=0, we find a luminosity 
distance of 3.5xl0^*cm, corresponding to an (> 20 keV) 
energy output of 1.5 xlO^^ erg for GRB 990510, based on 
the BATSE fluence of 2.56x10- " 
1999; see also Briggs et al. 2000). 

Unfortunately, the spectra taken on May 12 and May 16 
(not shown in Fig. 1) are of inferior quality, and cannot 
be used for absorption-line measurements. In the May 14 
spectrum, however, we detect Mgii again. The observed 
wavelength of Aobs = 7327 ± 5 A is consistent with the 
redshift determined in the first epoch spectrum. Its equiv- 
alent width is Wrest = 2.3 ± 0.6 A, which is (within the 
errors) identical to the measurement of May 11.2 1999 UT 
(0.8 days after the burst; Wrest = 2.6 ± 0.4 A). For the 
other absorption lines detected earlier it is not possible to 
obtain an accurate measurement of their EW. 

We will now derive a lower limit to the Hi column in 
the direction of GRB 990510 from the Fen lines at 2344 
A, 2600 A, and 2383 A in the first epoch spectrum. The 
relation between the EW of a line, W\, its column den- 
sity, Nj (number of atoms in the corresponding ionization 
state), and the oscillator strength, /, of the transition (j) 
is (Spitzer 1978): 

log ^ = \og{Nj X A X /) - 4.053 

Here the unit of W\ is A, the column density Nj is in 
cm^^ and the wavelength A is in cm. We do not use the 
strongest absorption line, Mgii A2800, since this line is eas- 
ily saturated in typical galaxy spectra, placing it on the 
fiat part of the curve of growth {W\ versus Nj x A x /). 
Even if the error in the EW determination were very small, 
this would still lead to a very uncertain value for the col- 
umn density. The Mgi line is not used either, due to the 
large uncertainty in the ratio of Mgi to Mgii, translating in 
a similar uncertainty in the column density. Fen, however, 
should be the dominant Fe ion in a dense neutral environ- 
ment. The following holds if the Fe lines are not saturated. 
Using the oscillator strength values from Morton (1991) for 
these lines (/(Fen A2344)=0.110, /(Fen A2383)=0.301 and 
/(Fen A2600)=0.224), and the values for Wa obtained from 
the spectrum of the first epoch (see Table ^) , we obtain 
log iV(Feii A2344) = 14.0 ± 0.1 cm'^, log 7V(Fen A2383) 
= 13.5 ± 0.2 cm-2 and log iV(Feii A2600) = 13.8 ± 0.1 
cm-2. We adopt the value log iV(Fe) = 13.8 ± 0.2 cm'^. 
Even though the EWs of the three different Fe lines result 
in similar values for the column, the lines are probably sat- 
urated, i.e. our estimate for the Fe column density should 
be considered as a lower limit. 

In converting the number of Fe atoms to a hydrogen col- 
umn density, we have to take into account that the metal- 
licity in high-redshift galaxies is likely to be lower than the 
Galactic value, and that a large fraction of the Fe atoms 
can be hidden in dust (Whittet 1992). A study of the 
metallicity in damped Lya systems (which have A^(Hi) ^ 



Table 2 
Absorption lines in GRB 990510 



Aobs (A) 



ID 



Zabs 



Wx. 



(A) 



4883 
5396 
6142 
6241 
6806 
7335 
7472 



Alni(1862.79)? 

Crll(2062.23)/Znll(2062.66)? 

Fell{2344.21) 

Feii{2382.76) 

Feii{2600.17) 

Mgii(2796.35/2803.53) 

Mgi(2852.96) 

weighted mean: 



1.6213 ± 0, 
1.6166 ± 0, 
1.6201 ± 0, 
1.6192 ± 0, 
1.6175 ± 0, 
1.6197 ± 0, 
1.6190 ± 0, 



0032 
0024 
0013 
0021 
0011 
0034 
0021 



1.6187 ± 0.0015 



0.8 ± 0.6 (0.010) 
0.8 ± 0.3 (0.010) 
0.6 ± 0.2 (0.011) 
0.5 ± 0.3 (0.012) 
0.9 ± 0.2 (0.011) 
2.6 ± 0.4 (0.017) 
0.6 ± 0.2 (0.013) 



Note. — W\ is given in the absorber's rest frame, i.e. the observed value is divided by (1+z). The error in W\ is dominated by the uncertainty 
in the placing of the continuum. The relative error expected from Poisson statistics and read noise is given in parentheses. The weighted mean 
redshift is calculated, using the separate measurements of all the lines, even though the identification of Aim and Cril/Znil is uncertain (see 
text). 
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Fig. 1. — VLT/FORSl spectra of the afterglow of GRB 990510; 0.8 and 3.9 days after the burst. The error spectra, calculated from the 
object and sky spectra using Poisson statistics, are also shown. Several absorption features (including telluric lines marked with ©) are 
indicated with the broken lines. 



102°cni"2) fj.Q^-^ ^ ^ 0.7 to 3.4 (Pettini et al. 1997; we 
adopt [Zn/H]=-0.8), allows us to estimate the Fe abun- 
dance at the redshift of GRB 990510 with respect to the 
solar abundance (Grevesse & Sauval 1999; log(Fe/H)0 — 
-4.5), obtaining log(iV(Fe)/iV(H))= -5.3. The Galactic 
Fe depletion factor in a cool disk environment is -2.2 dex 
(Sembach & Savage 1996, Savage & Sembach 1996), but 
here we adopt the typical depletion measured in damped 
Lya systems, which is around -0.6 dex. Taking all these 
points into account gives the following rough lower limit 
to the hydrogen column density: log Af(Hi) > 19.7 cm~^. 
Another and more robust way of estimating the Hi column 
density which is independent of dust corrections, is to use 
the Fell measured in DLAs (the gas phase only) around 
the redshift of the GRB. Using 10 systems with redshifts 
ranging from 1.2 to 2.0, we find [Fe/H] = -1.5 ± 0.5. As- 
suming that this is the most likely [Fe/H] abundance for 
the host galaxy of GRB 990510, we obtain log iV(Hi) > 
13.8 -I- 1.5 -I- 4.5 cm~^ = 19.8 cm~^, very close to our first 
estimate. A low column density is consistent with that 
found by Briggs et al. (2000) from fitting a combined set 
of BATSE and BeppoSAX X-ray and 7-ray data. 

HST imaging, performed in April 2000, appears to de- 
tect a very faint (V ~ 28) galaxy, located only 0.08" East 
of the position of the early afterglow (Fruchter et al. 2000, 
Bloom et al. 2000). This galaxy is probably responsible 
for the detected absorption lines in the spectra. The type 
and strength of the absorption lines, which are indicative 
of a low-ionization, high-density medium, strongly suggest 
that these originate in the host galaxy of GRB 990510. 
Even though the galaxy is extremely faint, the GRB af- 
terglow allows us to probe its interstellar medium and - 
in principle - to measure its metallicity. For the latter 
we need a measure of the column density that does not 
depend on the strength of the metal absorption lines (e.g. 
from a Balmer feature). 

4. THE EMISSION- AND ABSORPTION-LINE SPECTRUM OF 

GRB 990712 

Fig. 2 shows the spectra of GRB 990712, taken at 0.5, 
0.7 and 1.5 days after the burst (see Table n). The ob- 
vious emission lines are easily identified as [On] A3727, 
[Neiii] A3869, H7, HJ, and [Oiii] AA4959,5007. We also 
detect two absorption lines, which can be identified as 
Mgii A2796,2804 and Mgi A2853 at the same redshift as 
the emission lines, and are therefore intrinsic to the host 
galaxy. Using all these features in all three spectra, we 
obtain a weighted average redshift of z = 0.4331 ± 0.0004. 
In Table || we list the observed wavelength of the features, 
their identification and redshift, and the rest frame EW 
and flux (corrected for the Galactic foreground extinction), 
for the spectra taken at different epochs. The last two 
columns contain the EW and flux values averaged over all 
the spectra. 

It is important to investigate whether GRB host galaxies 
are indeed star forming galaxies, i.e. Hll galaxies, where 
massive O and B stars ionize the interstellar medium, giv- 
ing rise to prominent emission lines. These lines are also 
observed in galaxies that host an active galactic nucleus 
(AGN, e.g. Seyfert 2), where their strength does not only 
depend on the star formation, but on the nuclear activity 
as well. Most popular GRB progenitor models require a 
close connection with massive-star formation, but so far. 



this has not been conflrmed for any of the host galaxies. 

Hjorth et al. (2000) suggested that the host galaxy of 
GRB 990712 may be a Seyfert 2 galaxy on the basis of the 
[Oiii] A5007/H/3 ratio being greater than three (see Shuder 
& Osterbrock 1981), and its location in the log([Oiii] 
A5007/H/3) vs. log([Oii]/[Oiii] A5007) diagram (see Fig. 
2 of Baldwin, Phillips & Terlevich 1981). We also measure 
a ratio of [Qui] A5007/H/3 that is greater than three: 4.6 
± 1.0. However, from the more recent work of Rola, Ter- 
levich & Terlevich (1997), who employ the Canada-France 
Redshift Survey (CFRS) sample of emission-line galaxies 
at redshifts < z < 0.3, it is clear that a value of 4.6 is 
actually very typical for Hii galaxies (see their Fig. 1). 
Combined with our value for log([Oii]/'H/?) of 0.4 ±0.1, 
the host of GRB 990712 is located clearly within the Hll 
galaxy regime. We note that this diagram is corrected for 
extinction intrinsic to the distant galaxies, while the values 
in our Tableware not. However, the [Om] A5007/H/3 ratio 
is only slightly affected by reddening. Rola et al. (1997) 
also build non-extinction corrected diagrams to distinguish 
between Hii galaxies and other emission-line galaxies; in 
all but one of these, where its location is on the border 
of the Hii galaxy and Seyfert 2 regimes, the host galaxy 
of GRB 990712 is classifled as an Hii galaxy (e.g. Table 
I gives C3727-C4861 = -0.28 ± 0.03, whereas Fig. 4 of 
Rola et al. 1997 shows all definite Seyfert 2s have C3727- 
C4861 > 0.4), i.e. the emission lines arc produced by Hii 
regions that are being ionized by O and B stars. We con- 
clude GRB 990712 is most likely an Hii galaxy, and not a 
Seyfert 2. 

We have estimated the host galaxy extinction, by com- 
paring the observed ratio of H7/H/3 (0.26 ± 0.09) with the 
expected ratio for case B recombination (0.469 ± 0.009; 
Osterbrock 1989). Using the Galactic extinction curve of 
CardeUi et al. (1989), we obtain Ay = 3.4;^;7. The rela- 
tive large error is due to the marginal detection of H7. 

We can now estimate the star formation rate (SFR) in 
the host galaxy of GRB 990712 in three different ways: 
through the [On] and H/3 line luminosities (Kennicutt 
1998), and through the continuum fiux at 2800 A (Madau 
et al. 1998). For all these estimates a Salpeter initial mass 
function (IMF) has been assumed. Combining equation 
(3) of Kennicutt (1998), aluminosity distance of 6.8 x 10^^ 
cm (taking Hq = 70 km s"^ Mpc~\ f2n=0.3, and A=0), 
and the [On] A3727 line flux from Table g (the fluxes in this 
table have already been corrected for Galactic extinction) 
we obtain SFR([Oii]) = 2.7 ± 0.8 Mq yr-\ consistent 
with the value found by Hjorth et al. (2000). Taking into 
account the measured extinction (at Ha, due to the fact 
that the [On] SFR method is cahbrated through Ha), us- 
ing the Galactic extinction curve (CardeUi et al. 1989), we 
flnd SFR([Oii]) = 351:25^ M© yr'^ Using equation (2) of 
Kennicutt (1998), the case B (large optical depth) line ra- 
tio JHaJJHiS of 2.85 (Osterbrock 1989) and the H/3 flux of 
Table |, we find SFR(H/3) 1.7 ± 0.6 Mq yr"!. Corrected 
for reddening, this becomes SFR(H/3) == 64^^^° Mg yr"!. 
Finally, equation (2) of Madau et al. (1998), combined 
with the 2800 A flux (which is located at 4013 A at a red- 
shift of 0.4331), gives SFR(2800 A) = 2.8 ± 0.9 Mq yi-\ 
which becomes ~ 400 Mq yr^\ using Ay = 3.4. The [On] 
method, although it is indirectly calibrated through Ha 
and sensitive to the abundance and ionization state of the 



Table 3 
Absorption and emission lines in GRB 990712 









July 13.18 




July 13.42 




July 14.18 




Aobs (A) 


ID 


z 


W^A,rc.t (A) 


Line flux''''' 


Wx,r.st (A) 


Line flux 


Wx,rc.t (A) 


Line flux 


4015 ± 6 


Mgii AA2796,2804 


0.4340 ± 0.0021 


8.3 ± 1.4 




9.7 ± 1.9 




13.7 ± 4.5 




4096 ± 6 


Mgi A2853 


0.4357 ± 0.0021 


2.2 ± 0.6 




< 2 (2a) 




3.9 ± 1.4 




5342 ± 3 


On] A3727 


0.4332 ± 0.0008 


-47.1 ± 7.4 


3.27 ± 0.20 


-44.7 ± 7.9 


3.33 ± 0.22 


-45.8 ± 8.6 


3.37 ± 0.33 


5545 ± 5 


Neiii] A3869 


0.4332 ± 0.0013 


-8.4 ± 1.5 


0.60 ± 0.08 


-10.9 ± 2.3 


0.79 ± 0.11 


-7.5 ± 1.9 


0.55 ± 0.11 


6220 ± 3 


H7 A4340 


0.4332 ± 0.0007 


-8.1 ± 3.3 


0.50 ± 0.18 


-6.0 ± 2.5 


0.39 ± 0.14 


-4.5 ± 1.8 


0.29 ± 0.10 


6966 ± 3 


H/3 A4861 


0.4330 ± 0.0006 


-21.1 ± 3.9 


1.17 ± 0.12 


-28.6 ± 5.3 


1.59 ± 0.17 


-22.4 ± 7.7 


1.23 ± 0.31 


7106 ± 3 


Oiii A4959 


0.4330 ± 0.0006 


-43.7 ± 7.0 


2.36 ± 0.18 


-40.7 ± 6.2 


2.28 ± 0.16 


-42.8 ± 8.8 


2.38 ± 0.30 


7175 ± 2 


Oiii A5007 


0.4330 ± 0.0004 


-113.5 ± 17.1 


6.04 ± 0.31 


-110.7 ± 15.1 


6.15 ± 0.25 


-108.0 ± 18.6 


6.18 ± 0.53 




weighted mean: 


0.4331 ± 0.0004 















^All flux units: 10 ^® erg s ^ cm ^ 



''All line fluxes have been corrected for the Galactic foreground extinction. 
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Fig. 2. — VLT/FORSl spectra of the afterglow of GRB 990712, from 0.5 to 1.5 days after the burst. The error spectra, calculated from the 
object and sky spectra using Poisson statistics, are also shown. The absorption features (including telluric lines marked with ®) and emission 
features are indicated wit[i the broken lines. The observation dates (also shown in the figure), grisms, resolving power and signal-to-noise 
ratios are listed in Table 



""^^^(3 

I 



gas, is probably the least uncertain of the three. The H/3 
estimate is uncertain due to the unknown host galaxy stel- 
lar absorption underneath the emission. The SFR based 
on the ultra-violet continuum is very uncertain due to the 
fact that our flux calibration is an extrapolation below the 
V band, and, more importantly, the extinction correction 
at ultra-violet wavelengths is very uncertain. Normally, 
the 2800 A SFR is found to be lower by a factor of 2-3 (at 
redshifts up to z = 2.8) as compared to the Ha luminosity 
SFR method (Glazebrook et al. 1999, Yan et al. 1999). 

For the GRB host galaxies for which a SFR has been de- 
termined so far, even for the [On] method alone, the values 
range from 0.5 ± 0.15 Mq yr"! (GRB 970228; Djorgovski 
et al. 1999a) to 20 ± 9 Mq yr^^ (GRB 980703; Djorgov- 
ski et al. 1998); host galaxy extinction is not included, so 
these values should be considered as lower limits. Normal- 
ized by the host B band luminosity, this range is narrowed 
down to spread a factor of three. It has been noted (e.g. 
Djorgovski et al. 1998), that the range of SFRs for GRB 
host galaxies does not seem to be extra-ordinarily high, 
as compared to field galaxies at similar redshifts. E.g. 
Glazebrook et al. (1999) find a range of 20-60 Mq yr'^ 
(using Ha), for 13 field galaxies at z = 1, drawn from the 
CFRS. However, the galaxies at high redshifts for which 
these rates have been measured tend to be much brighter 
than the typical GRB host galaxy, and are therefore ex- 
pected to have much higher SFRs. The [On] equivalent 
width, which effectively is the star formation rate normal- 
ized by the blue band luminosity of the host galaxy, allows 
a more useful comparison. The mean Wrest ([On]) of the 
Glazebrook et al. sample is 33 with a standard deviation 
of 15 A, while we measure Wiost([Oii]) = 46 ± 1 A for 
the host of GRB 990712. This comparison also suggests 
that GRB 990712 occurred in a galaxy that is vigorously 
forming stars. 

Using the oxygen and H/3 emission lines, we can estimate 
the global oxygen abundance of the host of GRB 990712 
(see Kobulnicky, Kennicutt & Pizagno 1999, and refer- 
ences therein). From the values in Table we obtain 
-R23 = (^3727 + -^4959 + hmi) / H (3 = 8.9 ± 1.5, correspond- 
ing to log(0/H) = -3.7 ± 0.4 (we note that log(O/H)0 = 
-3.1; Cox et al. 1999). This estimate is based on fluxes not 
corrected for reddening, but such a correction would not 
change the abundance estimate substantially. For compar- 
ison: the oxygen abundances for a sample of 22 relatively 
nearby spiral galaxies range from log(0/H) = -2.7 to -3.5 
(Kobulnicky, Kennicutt & Pizagno 1999). 

5. THE TIME DEPENDENCE OF THE MGII FEATURE 

It is now well established that GRBs are the most en- 
ergetic events in the universe, with peak isotropic lumi- 
nosities up to 10^^ erg s~^ (Kulkarni et al. 1998). Also 
when the radiation from GRBs is beamed, the impact of 
the shock on the circumburst material is the same as in 
an isotropic explosion, and could be observable as time- 
resolved evolution of the Fe Ka line and edge in the X-ray 
regime (Weth et al. 1999, and references therein) and of 
ultra-violet (UV) absorption lines, redshifted to the op- 
tical domain (Perna & Loeb 1998). By monitoring the 
line evolution it is possible to obtain information on the 
density structure surrounding the explosion, and - in case 
the density can be measured by an independent method - 
the redshift to the burst may be obtained (Perna & Loeb 



1998). 

Our spectral observations of the afterglows of 
GRB 990510 and GRB 990712 extend over several nights, 
and both contain absorption lines, so that we can look for 
possible temporal evolution of these features. They are 
expected to decrease in time, if a considerable fraction of 
the atoms responsible for the absorption are in the vicinity 
of the site of the burst, and are ionized by the explosion 
(Perna & Loeb 1998). Alternatively, the burst may release 
atoms that are locked in the dust, which could result in a 
corresponding increase of the EWs. The clearest absorp- 
tion feature in both bursts is Mgii. In the May 11 and 
May 14 spectra of GRB 990510 (0.8 and 3.9 days after the 
burst, respectively), we measure Wiost — 2.6 ± 0.4 A (May 
11) and T4^rcst = 2.3 ± 0.6 A (May 14). For GRB 990712 we 
obtain Wrest = 8.3 ± 1.4 A, W^st = 9.7 ± 1.9 A and W^st 
= 13.7 ± 4.5 A for 0.5, 0.7 and 1.5 days after the burst, re- 
spectively. For both bursts, the values are constant within 
the errors. This is not very surprising; the Mgii feature is 
most likely saturated in the spectra of both bursts. This 
means that over a wide range of column densities, the EW 
is not expected to change by a detectable amount. On the 
other hand, a significant change in the EW would have 
indicated a large change in the column density. The Mgi 
feature in GRB 990712, which is possibly saturated as well, 
is also constant within the errors. 

6. CONCLUSIONS 

GRB afterglows allow us to probe galaxies that would 
otherwise be extremely difficult or impossible to study 
spectroscopically. We have determined a lower limit to 
the redshift of GRB 990510 through identification of sev- 
eral spectral absorption lines. The strength of the Mgi 
line is indicative of a cool, dense environment, which leads 
to the conclusion that the measured z = 1.6187 ± 0.0015 
most likely reflects absorption in the host galaxy ISM. 

Using both the absorption and emission lines in the 
spectrum of GRB 990712, we determine its redshift at z 
~ 0.4331 ± 0.0004. The emission-line ratios indicate that 
the host of GRB 990712 is an Hii galaxy, with an [On] star 
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formation rate (reddening corrected) of 35 
The large [On] equivalent width, compared to that of field 
galaxies at z < 1, also suggests that the host is vigorously 
forming stars. 

In order to put meaningful constraints on the circum- 
source medium, high resolution, high signal-to-noise spec- 
tra at several epochs after the burst are needed to resolve 
the velocity structures of the non-saturated lines, and al- 
low determination of the circumsource density distribu- 
tion and its evolution. This may become possible with the 
launch of HETE-II, since this satellite will allow follow-up 
observations of optical transients at early times, i.e. when 
they are bright (R ~ 16). Determination of the density 
profile could provide a major advance in solving the pro- 
genitor problem. 
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